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Abstract 

The hybrid meson is one of the most interesting new hadron specie 
beyond the naive quark model. It acquire a great attention both from 
the theoretical and experimental efforts. Many good candidates have 
been claimed to be observed, but there is no absolute confirmation 
about existence of hybrid mesons. 

In the present work we propose new calculations of the masses and 
decay widths of the hybrid mesons in the context of constituent gluon 
model. 

1 Introduction 

There is an intensive experimental activity attempting to detect new hadrons 
beyond the quark model which are predicted by QCD: glueballs, hybrids, 
diquonia, pentaquarks, ... 

Hybrid mesons (quark-antiquark-gluon) can have J^'" quantum num- 
bers which are not allowed by the naive quark model, like , O"*"", 1 '", 
2"'"~,...then they can't mix with the standard mesons and hence can facilitate 
their observation. 

These "exotic" objects are the most promising new species of hadrons 
allowed by QCD and subject of lot of works both in the theoretical and 
experimental levels. 
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From experimental efforts at IHEPt^l, KEK^l, CERNl^l and BNLW several 
J^'^=l^+cxotic resonances have been claimed to be identified, especially 
7ri(1600) and 7ri(1400) have received great interest, but some doubts are 
raised about the later one^^"^]. 

In the charm sector, the most interesting new states observed are the 
X(3872) and the Y(4260)[^l-. 

This work follows and completes several papers already devoted to this 
subject'^^^^l , new spin dependent estimations of the masses of hybrids and 
its consequences to the decay widths are proposed. 

In fact the hybrid mesons are studied from different models: lattice 
QCDl^^l, flux tube model'^^l, bag model'^^l, QCD sum rules'^^l and constituent 
gluon model'''"^^'^^"^^'. Some of them can perform both estimations of mass 
and decay widths. The nature of gluonic field inside hybrid is not yet be clear 
because the gluon plays a double role: it propagates the interaction between 
color sources and being itself colored it undergoes the interaction. Whereas, 
LQCD and Sum rules QCD make no assumptions about it, two important 
hypothesis can be retained from literature. The first one consider gluonic 
degrees of freedom as "excitations" of the "fiux tube" between quark and 
antiquark, which leads to the linear potential, that is familiar from quark 
model (flux-tube model). 

The second issue, which are supported by the present work, assumes 
that hybrid is a bound state of quark- antiquark and a constituent glue which 
interact through a phenomenological potential, precisely Coulomb plus linear 
potential supplemented by spin-spin , spin-orbit and tensor correction terms. 
The use of relativistic kinetics is appropriate for the study of the light flavor 
systems. The wave functions obtained are used for the calculations of the 
strong decay widths. 



2 Hybrid states and binding Hamiltonian 

For the representation of hybrid states we will use the notations: 

/g : is the relative orbital momentum of the gluon in the qq center of 
mass; 

Iqq : is the relative orbital momentum between q and q; 

Sqq : is the total quarks spin; 

jg : is the total gluon angular momentum; 

L '■ Iqq ® i g- 
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Considering the gluon moving in the framework of the qq pair, the Parity 
of the hybrid will be: 



P ^ P {qq) ■ P{g)- P (relative) 

= (-1) • (1) 

(—1) being the intrinsic parity of the gluon. Then the Parity of hybrid 
meson will be: 



P = + _ (2) 

The Charge Conjugation is given by: 

The QE and GE hybrids 

Sqq can takes values or 1; P and C impose Parity restrictions to Iqq, and 
Ig (Table 1). 

For lower values of orbital excitations ( Iqq and ^ 1) and 
P — -1, hybrid states can be built by two modes; Iqq — 1 and Ig — 0, 
which we shall refer as the quark-excited hybrid ( QEj, and Iqq = and Ig = 1, 
which we shall refer as the gluon-excited hybrid ( GE^. These particular modes 
play an important role in the decay selection rules to be mentioned later. The 
possible mixed QE-GE state is taked into account in this paper. 
The wavefunctions 

The total wave function can be written as: 

*S = ( ( (e., » ^€;) „ ® v^;;-') ® x2i) (4) 

^Jm(P,A) = {h^'^9^i^,\3,Mg)^lq^qq3gMg\Lm) 
{LmSqqIlqq \ J M) (5) 

The sum runs over the values of Iqq, mqq, Ig, nig, jg, Mg^ L, m, Sqg, jiqq 
and ^jLg excluding those not consistent with P and C. 



Here the Jacobi coordinates are introduced: 

^ 'a Mq + Mq ■ 

The Hamiltonian is constructed, containing a phenomenological potential 
which reproduces the QCD characteristics; its expression has the mathemat- 
ical "Coulomb + Linear" form, and we take into account also the additional 
spin effects. 

The basic hypothesis is the use of the relativistic Schrodinger-type wave 
equation: 

TV 



E + < + = E ^>{n). (6) 

In the case of multiparticle systems, another equivalent wave equation is 
more convenient l^^': 



where Mj are some "dynamical masses" satisfying the conditions: 
14// is the average over the color space of chromo-spatial potential: 



TV 



Veff = (v^)co;o. = ( - E • 

\ i<j=l I 

TV 

= E ^ij^i'^ij) ; (9) 
i<j=i 

where v{rij) is the phenomenological potential term. 
We take a potential which has the form: 

viuj) = vcLinj) + vsoinj); (lo) 

where the QCD-motivated " Coulomb. -I-Linear" term is of the form: 

vcL{rij) = -— + a Tij + c ] (11) 
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the as, cr, and c may be fitted by experimental data or taked from lattice 
and Reggc fits. 

The spin- dependent term can split into Spin-Spin, Spin-Orbit and Tensor 
terms : 

vsd^vss + vso + Vt (12) 



as f Si - X Pj Sj ■ X Pj 2 Si ■ Tij X - 2Sj • rj^- x pj 



(14) 

as f Si- Tij X Pi Sj ■ Tij X Pj 2 Si • Tij X - 2Sj • Vij x p 



i'"so)ij = 

(15) 

^^^^ij ^ ' ' - Si ■ Sj- ) ; (16) 

-'"i-'"j' ij 

The mass of the constituent gluon 

An important ingredient of the model is the mass of the constituent gluon 
irig. As for quarks, this represents a dynamical mass and is expected to add 
0.7 ~ 1.0 GeV to the corresponding quarkonia. The authors of ref. [17] 
are generating 800 MeV constituent gluon mass in the context of Dynamical 
Quark Model employing BCS vacuum, a value which is consistent with the 
(1600) glueball candidate. Furthermore, there is evidence for massive like 
dispersion relation for the gluon, with mass ranging from 700 ~ 1000 MeV , 
both lattice and from Schwinger-Dyson equations, f^^' 

The hybrid mass evaluation 

In earlier estimations'-'^'^l we have used nig — 800 MeV, here we adopt a 
slightly modified value (850 MeV) to the best fit of our results to the ones 
obtained by lattice calculations of l~^ccg and hbg masses. 

We have to solve the wave equation relative to the Hamiltonian: 



«=9> 9, 9 

with, for the hybrid meson : 



"9? 6 ' g 



(18) 
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The relative Hamiltonian is given by: 

pi A M„ mi Mg m? M„ mi 

with 

(20) 



+ Mq+Mq) ' 

The spatial trial wavefunctions are of the Gaussian-type: 



= exp Y;,(Q5) (21) 

For a J^*^ state we calculate the energy function: 

/«l»Jg(ftA).*S(p,A)dpdA 

and minimize it to respect of the parameters Mi, I3i^-i^ to have both masses 
and parameters of the corresponding wavefunctions. For more details see the 
Appendix. Numerical results are exposed in tables 2-4. 



3 The decay 

The decay process occurs through the quark pair creation mechanism, anni- 
hilating the constituent glue. Original quarks being spectators; more details 
can be found in [9] . The hybrid meson decaying in two standard mesons 
obeys to the following (model independent) selection rules: 

• only QE-hybrid meson decays in two S-standard mesons; 

• only GE-hybrid meson decays in one L and one S-standard mesons. 
Some earlier decay results are reviewed taking into account a significative 

mixing of the QE and GE modes generating from our mass estimations, this 
is summarized in tables 5-8. 
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4 Results and discussion 



The mass results are grouped in tables 2-4. 

We note that the hybrid masses obtained are in good agreement with 
those obtained by other models. The spin corrections are more important in 
the light sector, and do not exceed 19% for Aso and 7%. for Ass and Axens- 

The mass of the GE-hybrid is larger than those in the QE-modc. Indeed 
the strong force being proportional to the color charge, the exchange of a 
color octet does require an important energy. 

In the case of the pure GE mode, we have Mq-+ < Mi — < Mi-+ < Mq — 
(table 4). This is more close to the results of the "gluon excitation" models; 
this is simply because such works exclude the non excited glue {Ig = 0) mode 
(i.e. the QE one) from its construction of hybrid mesons. 

Note that l~'^nng is around 1.93 GeV (300 ~ 500 more heavier than 
7r(1600) and 7r(1400) candidates), this at least exclude that the last candidate 
to be hybrid meson'^1. 

On the other hand for = 1^ hybrids we find a significant mixing 
between QE and GE-modes. Using the l~~^nng wavefunctions obtained, we 
have estimate decay widths ( assuming Mi-+nng 1-6 GeV , Tables 5A-6A); 
the channel biir remains dominant, and the observation of I'^nng decaying 
in two S-mesons remains possible ( Vtot (1600)] ~ 60 MeV). However 
the mass gap imposes some doubts to consider the candidate 7r(1600) as an 
hybrid meson. This conclusion is supported ref. [6] . 

In the case for the candidate 1~+ around 2.0 GeV observed by AGS 
Coll. [^^1, which is more close to our mass estimations, (Tables 5B-6B), we 
note that for the "L+S" channel the decay widths are very large. We can 
observe a narrow resonance in the "S+S" channel. 

In the charm sector, the 1 ccg is estimated to have a mass around 
4.3 GeV which is consistent with the Y(4260) candidate. Otherwise, we 
find that 1 ccg decays in two (L-l-S)-mesons with decay widths sufficiently 

small to generate observable resonance ^^'^^ ~ 107 MeV^ . Note that 

the possible mixing of 1 ccg and the corresponding 1 cc is excluded '^1. 
Finally, to have an idea on the geometrical configuration of the hybrid 

mesons, we calculate the ratio |^ (table 9); the gluelump picture that con- 
sider the hybrid as a point like qq pair bound to a constituent gluon, is not 
supported by the present results. 
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Appendix: The matrix elements 

The calculation of the matrix elements requires the evaluation of many 
(cpu-intensive) six- dimensional integrals. Fortunately, the symmetry of the 
problem allows us to do this by an analytic way. 
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The spatial matrix elements 

Using the hyperspherical coordinates 

^ = (C; O5) = iC, fip; Qx; a) = (e; V^, 0x, 0) 
St'.; (0<^<f)and f = p^ + A^ (0<e<oo). 

the potential energy will be more symmetric and we can separate angular 
from hyperradial variables: 

VCp, A) ^ Vid) = + B^K] m ^ + c"; 

where [X] = (/p, rrip; Z^, ^a) = {Iqq, niqq; Ig, rrig) . 

Then matrix elements can be calculated analytically as: 

^ ^ 

where 

= C-J^'^e^'exp (-laf^j f ) ; = /p + /a 
The coefficients A[KiK2] and i3[i^^i^2] reads: 

with 

= n'^'' sin'^^cos'- ^ Yi^m,{^p)Yi^^,{nxy, 
dfl^ = (ifipdr^A sin^ 6* cos^ 6* (i^ ; dili = sin 6*2 (i6'j dcp^^. 

We have used the well known expressions: 

Z=0 m=-/ 

00 i 

Z=0 m=-i 
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where 



i+2 ; 
_ 1 1 rj, 

'^^ ~ 21+3;:^" 2^1^' 

r< = mm(ri,r2), 
r> = max(ri,r2). 

with f 1 = A and r2 = 

We calculate the integrals: 

oo ( 

z=o m=-; 
and we find 

''min 

X (ipiUlpi \y l(^-mpi-mp2) \ lp2fn) (^Al^Al | Yi(mAi+mA2) | ^A2^A2) , 

(/imi |Y,^| ^2^2) = jdn Yl^^{n)Yim{n)Yi,m,m; 

Imin = max(|/p2 - lpl\, \l\2 " ^Al|), 
^max = niin(/p2 + Ipl, l\2 + ^Al); 

where 

arctan(2) ^ 

/ sin2^cos2^rf^Vfe''^(^)^^^;^^(^)2r|g|r^+ / sin^^cos^^ 

axctan(2) 

Finally we can write 

2 
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In the same way we evaluate the S^^^^^j coefficient: 

2 





with 



("^ft )[/<'2i<ri] ~ <^(mpi+mp2,mAi+mA2)(~l) (T) 21 + 1 

I— I ■ 



where 

arctcLn(2) 



1 tan'+2^ 1 

■ + 



2/ + 3 2'+2 21-1 2^ 





2 



j sin^ ^ cos^ ^ d9 V^^f ^ (^) V^^;^- (^) (^|!^ cot'+2 9 - ^ cot' 9 



circtan(2) 



The terms A,^ ^ , and B,^ ^ , are expressed in terms of Gamma functions 
tions Bz{a,b) 



r (0), the "three J symbols" ( ) and the incomplete betha func 

' mi 7712 TTla ' 



The spin-spin parts 

In the same way we can calculate the spin-spin corrections using the 
expression: 



Svra/i ai ^ 1 1 , ^ . 
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where 



J^ 
J 



-I 



*jm(p, A)* expi-a'i p')^'^l^(p, A) dpdX, 



2 ^2\,TfPC, 



*5S(p,A)*exp(-a^, 
*jM(P,A)*exp(-a2^ 



X + x+p 



X — x-p 



)*;^(p,A) dpdX, 



We use here the developement: 



exp 



±pAcos(y, a")] = J](Tya(pA) P;[cos(y, a")] 



1=0 

oo 



4:71 



+1 



1=0 m=-l 



with 



+1 



2/ + 1 

Q(pA) = — - — / dp exp(pA p) Pi{p). 



Finaly we arrive to expressions containing a finite sums which are easy 
to evaluate. 

The spin orbit parts 

The spin orbit matrix elements can be expressed as: 
1 a. 



^Vso = 



6 2 



a L • L + a' 



-k) 
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with 



= 

A+ = 



c- 



Ml MnM-. 



Mg + Mg \Mg 
= -WgiWg^'^'^^^WS^' 

— — ^A"*" with q^q and g — > ; 

M, 2il^ 



Mg {Mg + Mg) \Mg 

(^B^ with q^q and g — >^ ; 



2 



X — + 



5, 



f / ' 



Mg {Mg+Mg 

^C^ with q ^ q and q ^ q^ ; 



M^ 



+ 



Mg{Mg + Mg)\{Mg + Mg) M, 



X ^ + 



9 / ' 



Mg {Mg+Mg 

D~ — — with q —>^q and q ^ q^ ■ 
We have to evaluate the following integrals: 



L 
K 



px^x 



1 



A ± x±p 



^^^(p,A) dpdX 
n X p'^.^'^^(p,A) (ipdA 



p or A. 
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As an example we calculate I^pp^'- 

1 



P X 



(-Va) ^^jg{p,X) dpdX 



A ± x±p 

The particular form of the trial wavefunction leads to expressions: 



I 



I± 



I 



x,y,z + circular permutation.; with: 



A ± x±p 



J *^^(P,A)* 



Pi 



A ± a;±p 
A. 



cxp ( -^"[W^]'^' ) *jm(p, A) dpdA, 



3 exp I ) ^jm(p> ^) oJpc^'^- 



A ± x±p 

These integrals can be calculated using the development: 

-(-3) _ 

~ 1 - .7:2 2 

=0 



y/l + x^± 2xp = Yl (t)' (2/ + 1) x'Pi (p); x^ ±1 

then we can write: 

te 



X± xp 



(-3) 



J- Jb 



1=0 l=-m 

If we have the same flavors, then Mq = Mg and 



te 

X 

In the end we have 



^ ^ < arctan 2 arctan 2 -< 6* -< | 
cos-3 e 2^ sin-3 e 



tan 9 



I± 



I± 



/ de sw? e cos^ 6 j^x^Att (-1)" {Ipi^pi InjY^^I /p2mp2) (/pimpi |njY;„| /p2"ip2) ; 



where n,- = 
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The tensor part 



The tensor corrections are calculated from the formula: 



6 M. 



+ 



+- 



[ ^ 3 {Sqi Sqj) Tij T {S q - S q)\ 



with the definitions 



PC, 



^■^;^(p,A) dpdX, 



A ± a;±p 
1 



n± = 



A ± x±p 
A ± x±p 



*jm(^,A) dpdX, 



A ± x±p 



To evaluate these integrals we need to exploit the development: 



A ± xp 



(-5) 



J2 (t)' [21 + ?,- {21 - 1) x'] 47r YiMYimi^x), 



l=—m 
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p 


c 




Iqq 


h 









even 


odd 






1 


odd 


even 




+ 





odd 


even 




1 


even 


odd 









even 


even 


+ 




1 


odd 


odd 









odd 


odd 


+ 


+ 


1 


even 


even 



Table 1: Iqq and Ig Parity. 





1 — 


0-- 


0-+ 


1-+ 


nng 


1.40 


1.58 


1.73 


1.93 


ssg 


1.66 


1.87 


2.02 


2.21 


ccg 


4.27 


4.35 


4.38 


4.48 


bbg 


10.50 


10.66 


10.68 


10.80 



7116/6 2: Mixed QE-GE hybrid mesons masses with spin corrections (in GeV) 





QE 


GE 


nng 


1.31 


1.70 


ssg 


1.57 


2.00 


ccg 


4.09 


4.45 


bbg 


10.34 


10.81 



Table 3: Pure QE and GE-hybrid masses (in GeV) 





0-+ 


1 — 


1-+ 


0— 


nng 


1.83 


1.90 


2.04 


2.06 


ssg 


2.13 


2.21 


2.33 


2.38 



Table 4- Pure GE-light hybrid masses with spin corrections (in GeV) 



17 



r 


1.5 


r 


53 




7.0 


Ttot (1600)] 


61 



Table 5A: Decay widths of the (M=1.6) hybrid in (S+S)-standard mesons (in agMeV). 





31 


pn 


100 


rp(1450)7r 


18 


Tk*k 


32 


rA'*(1410)X 


0.3 


Tu.L (2000)] 


181 



Table 5B: Decay widths of the (M=2.0) hybrid in (S+S)-standard mesons (in asMeV). 



L 





1 


2 




145 


433 


722 




145 


434.7 


724.5 


r/0(1285)7r- 


117 


88.2 


147.4 


r/p(l420)7r- 


39 


29 


47.9 


rtoi[i-+ (1600)1 


446 


563 


1642 



Table 6A: Decay widths of the (M=1.6) hybrid in (L+S)-standard mesons (in asMeV). 



L 





1 


2 




477 


1425 


2378 


r/0(l285)7r- 


481 


361 


604 


r/0(1420)7r- 


524 


340 


643 



Table 6B: Some partial decay widths of the (M=2.0) hybrid in (L+S)-standard mesons (in agMeV). 
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L 





1 


2 




129.5 


388.5 


647.5 


Td+d- 


135.1 


406 


676.2 


^ D+D7 


142.8 


428.4 


714 


r — r 


0.00 


0.00 


0.00 


s = o 

^D*+D*-=^]:)*OD*o "S" = 1 

5 = 2 


30.8 
0.00 
49.3 


92.4 
0.00 
369.6 


1.47 
0.00 
24.5 


Ttot [I— (4.26)] 


982.1 


2463.3 


2914.1 



Table 7; Decay widths of the (M=4-26) hybrid in (S+S)-standard mesons (in agMeV). 



£)i(2420)D" ~ -^r'i(2420)D" 
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Table 8: Partial decay widths of the (M=4.3) hybrid in (L+S)-standard mesons (in agMeV). 





1-+ 


0-+ 


1^^ 


light sector 


0.8 


1.0 


1.4 


h-sector 


0.6 


0.8 


1.1 



Table 9: Geometrical configuration of the hybrid mesons. 
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